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Abstract 

Rapid advances are being made toward optically cooling a single mode of a micro-mechanical sys- 
tem to its quantum ground state and observing quantum behavior at macroscopic scales. Reaching 
this regime in room-temperature environments requires a stringent condition on the mechanical 
quality factor Qm and frequency Qmfm^kBTha.th/h, which so far has been marginally satisfied 
only in a small number of systems. Here we propose and analyze a new class of systems that should 
enable unprecedented Q-frequency products. The technique is based upon using optical forces to 
"trap" and stiffen the motion of a tethered mechanical structure, thereby freeing the resultant 
mechanical frequencies and decoherence rates from underlying material properties. 



The coupling of a high-Q mode of a micro-mechanical oscillator to an optical cavity has 
emerged as a promising route toward observing quantum behavior at macroscopic scales [T] . 
This opto-mechanical interaction is being used, for example, to optically cool a mechanical 
mode toward its quantum ground state [Ij. Ground-state cooling requires that the product 
of the quality factor Qm and frequency fm of the mechanical mode exceed ksTxyath/h, where 
h is Planck's constant. For a room-temperature bath, this condition is marginally satisfied 
only in a small number of current experiments [2]. The ratio /iQm/m/^B^bath also determines 
the quantum coherence time of the system relative to the mechanical period. Significant 
improvements to Q-frequency products are thus critical for schemes to prepare and detect 
non-classical states of motion |3], H] . Mechanical systems exhibiting extremely high quality 
factors also offer novel opportunities for precision measurement and force detection |5],|6]. 

In this Letter, we propose a class of systems that should enable unprecedented Q- 
frequency products. The approach is based upon optically "trapping" a tethered membrane 
with low natural mechanical frequency in the anti-node of a strong optical standing wave (see 
Fig. [l^) ^ . While there are many possible realizations, here we focus on a pendulum ge- 
ometry, where a relatively large disk is supported by a single thin tether. The dielectric 
disk is attracted to the anti-node of the field, leading to an optical stiffening of its flexural 
modes. Under realistic conditions, the re-normalized mode frequencies can be significantly 
enhanced over the values expected under material stresses alone. Of particular interest is 
the "center-of-mass" (CM) mode, where the disk oscillates in the optical potential with 
negligible flexural motion. We show that this motion exhibits an extremely large ratio of 
potential energy stored in the optical field to strain energy, f/opt/t^mech- This can yield a 
correspondingly large increase in the Q-frequency product over a conventional mechanical 
system due to the suppression of dissipation through internal friction. 

Our approach to achieving long coherence times builds upon previous proposals, which 
suggested that the highly isolated CM mode of an optically levitated nanosphere P [H] can 
enable quantum opto- mechanics in room-temperature environments [inilll]. Compared to 
the nanosphere, our approach has two significant advantages. First, the nanosphere scat- 
ters light omni-directionally, leading to motional heating via photon recoil. Suppression of 
recoil heating to reach the quantum regime requires spheres with sub-wavelength volumes, 
[TU]. In contrast, the planar membrane primarily couples the counter-propagating 
components of the trapping beam, strongly reducing recoil heating even for large systems. 
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Second, these membranes can be fabricated using well-established techniques that already 
yield excellent mechanical and optical properties in a number of experiments [21 [12] • Our 
proposal thus shows how the ideas of optical levitation can be brought to bear upon "conven- 
tional" and practically deployable mechanical systems to yield remarkable coherence times. 

We begin by considering the mechanical modes of a free, thin circular membrane with 
thickness d and radius a in the absence of any tethers. In equilibrium, the membrane is 
situated at 2; = 0, in the anti-node of an optical standing wave E(z) oc cos kz (see Fig. [T^). 
The optical field polarizes the dielectric disk, yielding a gradient force trap around z = 0. 
Absent any internal forces, the optical field traps a thin membrane [d ^ A) with a restoring 
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frequency given by uJopt{r) = ^^l!c'^~^^ ^ ' where r is the radial coordinate. Here /(r) is 
the beam intensity profile (assumed to be rotationally symmetric) in the direction transverse 
to z, = 27r/A is the optical wavevector, p is the mass density, and e is the dielectric constant. 
Now including the internal stresses, the mechanical displacement field from equilibrium, 
C{x,y), obeys 

subject to free boundary conditions [13j (also see Appendix). E,a denote the Young's 
modulus and Poisson's ratio, respectively. Here and in the following, is understood 
to be the Laplacian in the transverse plane. Due to the rotational symmetry, we seek 
solutions of the form ({x,y) = /(r) cos m6'e~*'^'"*. The spatial modes are indexed by the 
number of nodal diameters and circles, (m, n) (see Fig. [I]d). For our numerical results we 
take material parameters E = 270 GPa, a = 0.25, p = 2.7 g/cm^, e = 4 corresponding 
to stoichiometric silicon nitride [2J, and an operating wavelength of A = 1 pm. For a 
free disk without optical forces, the fundamental (2, 0) fiexural mode has natural frequency 

(2,0) /o^^n oc; d I E _ , ,(2,0) 



^rn,'nat/2^^0-25^y p^^_^2) , or uj)^'J^^/2ti^1.?) MHz for a disk of dimensions a = 10 /im, 
(i = 50 nm. There is also a trivial solution corresponding to the CM or (0, 0) mode with 
constant /(r) and zero frequency. 

For a uniform intensity, /(r) = Jo, the natural radial functions f'i!]^t^\r) remain eigen- 



modes of Eq. (A3), but with re-normalized mechanical frequencies given by ujln''"\l{ 
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\J uj1p^^{Iq) + (w^nat)^- Thus, the CM mode is now a non-trivial solution with frequency 
= Uopt, but still retains a uniform spatial profile. The frequencies of all the fiexural 
modes increase as well, with the CM mode remaining the lowest in frequency (see Fig. [2| 
The absence of energy stored in internal strains for the CM motion has the important 



implication of eliminating dissipation due to internal friction. Instead, the energy is stored 
in an optical potential that contributes no losses (but can contribute a recoil heating force, as 
described later). To quantify this, we consider the influence of thermoelastic damping on our 
system. We focus on thermoelastic damping because a) it can be analytically modelled [131 
S], b) it is a fundamental limit even for "perfect" devices [H], and c) a number of micro- 
mechanical systems are approaching this limit [151 US]- We emphasize, however, that our 
conclusions are qualitatively correct for any internal dissipative process. 

Thermoelastic damping arises because realistic materials have a non-zero coefficient of 
thermal expansion. The flexural motion creates local volume changes that then lead to 
temperature gradients and heat flow. Mechanical energy must be expended to drive this 
heat flow, leading to a flnite Qm- We make two simplifying assumptions to the general 
thermoelastic equations [TSl [T^ . which are well-justifled in our system. First, we assume 



that the thermoelastic coupling is weak, so that the spatial modes determined by Eq. (A3) 
are not altered to lowest order. Second, the strains vary most rapidly along the thin direction 
of the disk, and thus we ignore the relatively small transverse temperature gradients. The 
temperature fleld is given by T{x,y, z,t) = Tbath + ^T{x,y, z,t), where AT satisfies the 
driven heat equation 

with boundary conditions d/ST/dz = at 2; = ±(i/2. Here cy is the heat capacity per 
unit volume, n^h is the thermal conductivity, and a is the volumetric thermal expansion 
coefficient (we take cy = 2 J/cm^-K, = 20 W/m-K, a = 4.8x10^*^ K^^ for SiN). The 
work done in driving the heat flow over one cycle is 

AW^--^ / dt dhAT{r){d^AT/dz^), (3) 

bath JO J 

and the thermoelastically limited quality factor is Qm,th = 27r(f/opt + Uraech)/AW, where 
f/opt and t/mech are the energies stored in the optical field and strains, respectively. To good 
approximation, one finds that the thermoelastically limited Q-frequency product is given by 

Qm,thfm 

7rEd22Tl''tha2 i+f (-'- + ^opt/f^moch) (scc Appendix). Thus the storage of energy in the 
optical field leads directly to an enhancement of the Q-frequency product. As a useful com- 
parison, in the absence of optical trapping, an unstressed stoichiometric SiN film of thickness 
d = 50 nm would have a Q-frequency product limited to Qm,th/m~4x 10^^ Hz at room tem- 
perature, which only marginally exceeds the fundamental limit /csTbath/^^Sx 10^^ Hz needed 



for ground-state cooling. The Q-frequency product is also proportional to the number of 
coherent oscillations that the system can undergo before a single phonon is exchanged with 
the thermal bath, A''^.''^'^'^^ = Qm,thfmh/{27rkBTi,a.th) i^l^^'^^^l for the conventional membrane 
described above). A large value is critical to preparing and observing quantum superposition 
or entangled states |3], H]. 

We now examine thermoelastic damping of our free disk. Clearly, if Wopt(^) is spatially 
uniform, the CM mode has no internal strains and experiences zero thermoelastic dissipa- 
tion. A finite beam waist creates inhomogeneous optical forces that internal stresses must 
compensate for, which mixes CM and internal motion together (we still refer to this mixed 
mode as the "CM"). The mixing becomes significant when the variation in Uopt{r) across the 
disk overtakes the natural rigidity of the system (as characterized by the natural fundamen- 
tal frequency), and can be avoided by using sufficiently large beam waists or low intensities. 
In this regime, significant enhancements to the Q-frequency product should result. 

These effects are illustrated in Fig. |3^. For concreteness, we assume that the trapping 
beam has a Gaussian profile with waist w, I{r) = Iqc'"^^^^^^ (for now, we ignore possible 
corrections due to distortion as the beam diffracts around the disk). In Fig. [3^, we plot the 
ratio f/opt/f^iech for the CM mode as a function of its frequency, which is varied through the 
peak intensity Jq. We use the same disk dimensions as before and a waist of w = 35 fim. The 
energy ratio monotonically decreases, reflecting the increased inhomogeneity in the optical 
potential. In the inset, a complementary process is illustrated, where the beam waist w is 
varied, while the peak intensity Jq is flxed such that a;opt(/o)/27ir = 1 MHz. The energy ratio 
increases indeflnitely with w/a and approaches inflnity in the plane- wave limit. 

We now consider the realistic pendulum geometry shown in Fig. [l^, where the tether 
provides an extremely weak restoring force for the "CM" motion of the disk. We flrst 
present a simplifled analysis that isolates the role of the tether on the mode spectrum and 
Q-frequency product. Speciflcally, we treat the membrane as a perfectly rigid point particle 
of mass M, which experiences an optical restoring force with frequency Uopt, while internal 
stresses alone act on the tether. Then, for a tether of length L whose long axis is situated 
along X, the displacement fleld (f){x,t) (where < x < L) satisfles the beam equation [13], 

Here h denotes the width of the tether (assumed to be square in cross-section). The beam is 
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clamped at x = 0, 4>{0,t) = dx4>{0,t) = 0, while at x = L the boundary conditions are given 
by dl(j){L,t) = and Mdf(j){L,t) = -Mujlp,(f){L,t) + Eb^dl(f){L,t)/ 12. The last equation 
describes the acceleration of the membrane due to optical restoring forces and the shear 
force imparted by the tether. 

It is straightforward to solve for the system eigenmodes (see Appendix) and the re- 
sults are summarized here. For large mass ratios between the membrane and tether, 
M/mj— >oo, the modes usually consist of a CM mode for the membrane with frequency 
~ + Wopt and a set of discrete tether modes with frequencies u;„ ~ (^ ^"^^^^^ j > 
where /3 = {Vlp/ ElP'Y^^. The CM mode spectrum is understood as a low-frequency "pendu- 
lum" mode (with natural frequency ujp ~ a/ Eb'^/AML^, where cUp^Wopt, for our systems 
of interest) whose frequency can be strongly re-normalized by the optical force, while the 
tether mode spectrum results from the heavy membrane essentially acting as a second clamp. 
This description holds except near degeneracies Uopt^oJn, where coupling between the tether 
and membrane motions yields an avoided crossing whose width decreases with increasing 
mass ratio M/m^. This result is illustrated in Fig. [2|d, for a mass ratio of M/rrit = 125 (cor- 
responding to the disk size considered earlier attached to a tether of length L = 50 fim and 
width b = 50 nm). In Fig. |3}d, we plot the energy ratios f/opt/f^mech for the CM motion as 
a function of Uopt- Here, the strain energy is completely attributable to the tether, as we 
take the membrane to be a rigid object. The energy ratio is dramatically reduced near the 
avoided crossings, while away from these crossings, the energy ratio plateaus to a value near 
Uopt/Uraech^^M/mt (see Appendix). 

For a realistic tethered system (as in Fig. [i]d) where the membrane is not perfectly rigid, 
mode mixing between the tether and membrane and mixing between the CM and internal 
membrane motion will occur simultaneously. We have numerically solved the full stress- 
strain equations for such a system using COMSOL, a commercial finite-element simulation 
package. A characteristic mode spectrum is plotted in Fig. [2]:3 as a function of the peak 
trapping intensity Jq, for parameters a = 10 yum, L = 50 fim, d = b = 50 nm, and 
w = 35 /im. Away from avoided crossings, the modes can clearly be characterized as 
tether modes (gray points) or membrane modes (color). For our particular choice of beam 
waist size, the tethers themselves experience optical restoring forces, leading to a slight 
optical stiffening of tether modes that have displacements along the optical propagation 
axis. Comparing the membrane modes, the CM mode lies lower in frequency than the 
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flexural modes, as in the case of a free disk. A nearly degenerate torsional (1,0) mode also 
exists, which in principle should have no opto-mechanical coupling to the cavity field and 
can be ignored (see Appendix). 

In Fig. we plot the energy ratio t/opt / f^mcch for the CM mode of our tethered structure 
as a function of its frequency. The features displayed here are clearly a combination of 
those appearing in the limiting cases of a free disk and a rigid disk attached to a tether. 
In particular, large plateaus in f/opt/f^mcch appear away from avoided crossings with tether 
modes. The plateau heights are higher than the limit f/opt/^mcch~8M/mt predicted by the 
simple model, since the tether experiences an optical spring force as well. The decrease in the 
plateau heights with increasing CM frequency is associated with increased mixing between 
pure CM and internal membrane motion. For the realistic geometry considered here, an 
enhancement in the Q-frequency product on the order of ~10^ compared to a conventional 
system can be realized at a frequency of Co'm/2vr~l MHz. For a thermoelastically limited 
system, this corresponds to a coherence time of iV^^^'^'^''~10^. 

Thus far, we have assumed that the trapping beam has a Gaussian profile. For cavity 
opto- mechanics p], it will be necessary to trap the membrane within a Fabry-Perot cavity, 
as illustrated in Fig. [l^. For example, here, a relatively strong beam could be used for trap- 
ping, while a second, weaker beam with a non-zero intensity gradient at the trap position 
would facilitate cooling of the CM motion or quantum state transfer processes [TT]. The 
membrane scatters and diffracts the cavity light, which introduces two important effects. 
First, the mode will no longer be Gaussian, and the new optical mode accommodated by 
the cavity mirrors and the corresponding optical forces must be determined. Second, photon 
scattering out of the cavity reduces cavity finesse, and the associated random momentum 
kicks ( "photon recoil" ) imparted on the membrane lead to additional decoherence. We have 
performed detailed simulations of the modified cavity fields and coupled them to the equa- 



tion of motion (A3) for the membrane (see Appendix). We find that for realistic cavity 
geometries, a value of N^^^^^^IO"^ in a room-temperature environment can still be obtained, 
taking into account both thermoelastic processes and recoil heating. Remarkably, the co- 
herence time of this system is comparable to that of a much smaller levitated nanosphere of 
radius r~25 nm [10], or a conventional 1 MHz oscillator with Qm exceeding 10^ at a bath 
temperature of 1 K. 

We emphasize that our approach to reaching the quantum regime is fundamentally dif- 
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ferent than other "optical spring" proposals based upon optical backaction forces [T71 [18] . 
In the latter case, the linear coupling of the mechanical displacement to the intra-cavity 
intensity can yield a dynamic optical spring effect. This effect, however, is accompanied by 
significant Raman scattering of the optical pump field, which causes phonons to be rapidly 
removed and added to the system. While this does not preclude ground state cooling [18], a 
more detailed analysis (see Appendix) shows that it imposes severe limitations on the quan- 
tum coherence time and makes it difficult to prepare, e.g.^ quantum superposition states. 
In contrast, in our scheme, the phonons are truly long-lived excitations. Furthermore, re- 
gardless of the trapping scheme used, our analysis properly captures the role that strong, 
spatially non-uniform optical forces have in mixing internal motion, which is neglected in 
lowest-order opto-mechanical models but relevant to most flexural systems (see Appendix). 

Although we have focused on thermoelastic losses in the above calculations, we expect 
similar improvements for any other internal damping mechanism. The key idea is that it 
is possible to circumvent natural material limits of damping by storing energy in a lossless 
optical field rather than the internal strain. By making the ratio of these energies large, 
f4pt/f^mcch^l, any internal losses can be suppressed by a corresponding degree. This fun- 
damental observation allows one to design a novel class of mechanical systems that can 
be fabricated and deployed using conventional techniques, yet yield Q-frequency products 
that are several orders of magnitude higher than previous systems. We believe that this 
work will stimulate further investigation into the relationship between optical forces and 
material dissipation in a number of systems where the mechanical motion can be strongly 
renormalized by light [18], HQ]. Furthermore, we anticipate that such studies will open up in- 
teresting possibilities for quantum manipulation of mechanical systems in room-temperature 
environments. 
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Appendix A: Equation of motion for a free disk 



Here we derive the equation of motion for a thin non-uniform disk of thickness d{r), 
which has a reflection symmetry around z = (such that the surface of the plate is located 
at z = ±(i(r)/2). We are interested in the situation where the thickness is much less than 
the characteristic transverse size {e.g., the radius a of a circular disk), such that its degree 
of freedom along the thin direction can be effectively eliminated and the flexural motion 
can be described by a two-dimensional displacement field ({x,y). The equation of motion 
for ({x,y) can be obtained by a generalization of the derivation for a uniform disk given in 
Ref. [13] • Specifically, the energy associated with the displacement field ({x,y) is given by 

E f 

t^mech = ^^^^ _ J dxdyd{x,yf [{C. + Cyyf + 2(1 - a) {{C^yf - CxxCyy)] ■ (Al) 

Here E, a are the Young's modulus and Poisson's ratio, respectively, while C,xx = ^i etc. To 
derive the equilibrium field C,{x, y) under some external normal pressure P(x, y), we employ 
the variational principle to minimize the system energy. Under small variations 5^, and 
following some algebra, the variation in f/mech can be written as the sum of an integral over 
the transverse area of the disk and two integrals over the circumference or edge of the disk, 

5f/mech = J dxdy Zi5C + j dl + j dl Z^dC- (A2) 

Here Zi are complicated expressions involving C, and d{r) whose forms are given below, while 
n denotes the normal to the edge of the disk. The integral over the disk area yields the 
equilibrium equation of the disk, Zi = P{x,y), or the dynamical equation can be obtained 
by replacing P{x,y)-^ — pd{r)^. Doing so, and including the effect of external optical 
forces, one finds 

(A3) 

Here we have defined g{r) = d{r)^, and is understood to be the Laplacian in the trans- 

1 /2 

verse plane. As described in the main text, (X'opt(?") = • a disk with free 

boundary conditions at the edge, the quantities 6( and d5C,/dn are arbitrary on the bound- 
ary, so the coefficients ^2,3 should vanish, yielding the two boundary conditions. Defining n 
and I to be the normal and tangential directions to the edge of the disk, and 6 as the local 
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angle between x and n, these boundary conditions become 

d 



= ^ + {l-a) 



\ gCxy COS 26 + {g/2) sin 26'(C^^ - (yy)) 



dn 

+ cos 6{Cyygx - CxyGy) + SmO{Cxxgy " CxyQx)] , 

= V^C + (1 - ^) [Kxy sin cos - Cyy cos^ 9 - (^^ sin^ 6) . (A4) 
Appendix B: Thermoelastic damping 

We begin with Eqs. (2) and (3) in the main text that describe thermoelastic damping, 



and 

p2lT/oJm r 

AW^-— / dt dhAT{r){d^AT/dz^). (B2) 

^bath Jo J 

The first equation describes the driven heat equation along the thin direction of the mem- 
brane, with boundary conditions dAT/dz = at 2; = ±d{r)/2. Here cy is the heat capacity 
per unit volume, Kth is the thermal conductivity, and a is the volumetric thermal expansion 
coefficient (we take cy = 2 J/cm^-K, = 20 W/m-K, a = 4.8x10"^ for SiN). The 
second equation describes the amount of work done in driving the heat flow over one cycle 



for a particular mechanical eigenmode. Solving Eq. (Bl) and substituting into Eq. (B2) 
one finds 



losoMi-^)^ J ^ ^ 

Let us now compare this quantity with the total strain energy f/mech given in Eq. (Al). 
For simplicity, here we specialize to the case where the disk has a uniform thickness d, such 
that 

Ed^ r 

t^mech = l^^^^—T^ J dxdy [iV'Cr + 2(1 - a) {{CxyY - CxxCyy)] • (B4) 

Note that AW and f/mech have similar forms, as both involve an integral over the membrane 
area of the quantity (V^C)^- The strain energy contains a second term, however, whose 
relative importance we characterize now. The second term can in fact be re-written as a 
line integral around the circumference of the membrane, 



/ 



dxdy{Cxyf - CxxCyy = f dx CxyCx " f dy CyyCx- (B5) 
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This boundary integral vanishes identically for certain types of shapes or boundary con- 
ditions, such as a clamped membrane. For our free disk, this boundary term does not 
identically vanish, but numerically we can confirm that the boundary contribution is 
small relative to the total strain energy. To good approximation then, we can write 
f^nech ~ 24(f-g-2) / '^^ '^V (^^C)^- This Icads to the expression for the thermoelastically lim- 
ited Q-frequency product given in the main text. 



Qm,thfm - Tpiorr 2l~, ^ + 77 " (^^) 



Appendix C: Rigid membrane attached to tether 

Here, we derive in detail the properties of a rigid membrane of mass M attached to a 
single tether. As described in the main text, the tether has a length L (oriented along the x- 
axis) and a square cross-section of width b. The displacement field t) (where < x < L) 
satisfies the beam equation | il3j . 

The beam is clamped at a; = 0, 0(0, t) = 9a;0(O,t) = 0, while at x = L the boundary 
conditions are given by dxx4'{L, t) = and the force equation of the membrane attached 
there, Md^(f){L,t) = —Mujl^^(j){L,t) + Eb^dxxx4>{L,t)/ 12. The last equation describes the 
acceleration of the membrane due to optical restoring forces and the shear force imparted 
by the tether. 



The boundary conditions at a; = require that the general solutions of Eq. (CI ) take the 
form 

0(x) = ci(sin kz — sinh kz) + C2(cos kz — cosh kz), (C2) 

where the dispersion relation is given by w = {k/f3Y and (3 = {12p/ Eb'^Y^^. The two 
boundary conditions at x = L can be written in matrix form as Q{ci €2)^ = 0, and the 
corresponding equation for the eigenfrequencies, det Q = 0, reads 

M{uj'^ — J (cos 7 sinh 7 — sin 7 cosh 7) H — (1 + cos 7 cosh 7) = 0, (C3) 

where 7 = fiL^/Zj. With no optical forces (wopt = 0) and large membrane to tether mass ratio 
Mjm{^\^ the solutions to the above mode equation consist of a low-frequency "pendulum 
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mode" with frequency oOp ~ ■>/ Eb'^/AML^ and a set of discrete tether oscillation modes 
with frequencies Un, which approximately satisfy the relation cos7sinh7 — sin 7 cosh 7 = 0. 
For sufficiently large n {i.e., large enough 7), the natural tether frequencies asymptotically 



approach Un = (7m + 7r/4) . With optical forces, a power series expansion of Eq. (C3) 



reveals that the pendulum mode frequency becomes re-normalized to the value u. 



CM 
m 



Up + Wopt, which can be associated with the CM mode. 

Generically, for sufficiently large optical forces Uopt S> Up, a large membrane mass 



M means that the eigenvalue condition of Eq. (C3) is approximately satisfied when 
M(co'^ — (X'qp^.)(cos7 sinh7 — sin 7 cosh 7) ^ 0, which means that the solutions typi- 
cally consist of a single CM mode with frequency Um^uJopt or tether modes satisfying 
cos7sinh7 — sin 7 cosh 7 ^ 0. When both conditions are satisfied simultaneously {e.g., 



near a degeneracy point), the second term of Eq. (C3) must be taken into account, which 
yields avoided crossings that mix the tether and CM modes together. 

We now consider the ratio of energy stored in the optical field to the strain energy, 
f4pt/f^mcch, for the CM motion away from an avoided crossing. The energy stored in the 
optical field for this system is simply given by f/opt = (l/2)^i^opt</'(-^)^5 while the strain 
energy in the beam is given by [13j 

It can readily be shown that for frequencies where 7^1, the spatial modes are given to good 
approximation by 

(f){x)^ci (sin kx — cos kx + e"*""") . (C5) 

Then, when the CM mode is positioned halfway in between two tether modes (say 
at a frequency u^^{7in + 37r/4)^/(/3L)^), evaluation of the stored energies yields 



Appendix D: Dynamic optical spring 

An increase in the frequency of a mechanical mode can also be achieved through a dy- 
namic backaction effect in an opto-mechanical system [T71 HE], as opposed to the "static" 
optical potential considered in our work. We briefly describe the dynamic effect and com- 
pare the two mechanisms here. Speciflcally, we consider a mechanical degree of freedom 
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whose displacement is linearly coupled to the optical cavity frequency. The corresponding 
Hamiltonian for such a system in a rotating frame is fIU[ fIT\ 

Hint = 1 — i^^m^"^ ~ ^Sa'd — hhj'za^a — h^lL{d + d^). (Dl) 

2m 2 

Here d is the annihilation operator for the optical mode, z,p are the position and momen- 
tum operators corresponding to the mechanical resonator, Um is the natural mechanical 
frequency, 6 = ul ~ ouq is the frequency detuning between an external pump field driving 
the optical cavity and the optical resonance frequency uo (when the mechanical resonator 
is in equilibrium), Ql is the driving amplitude, and u' is the optical cavity frequency shift 
per unit mechanical displacement. In addition to the Hamiltonian terms, the optical cavity 
is assumed to have losses characterized by a linewidth k. 

For weak opto-mechanical coupling, it is customary to linearize the optical cavity dy- 
namics around the classical steady-state value (d) = a = iVLL/{n/2 — i5) (here we have 
incorporated a steady-state shift of the optical resonance frequency into our definition of 
the detuning 5), and eliminate the cavity to yield an effective susceptibility xi!^) of the 
mechanical displacement in response to an external force /(w) |2D]- Specifically, one finds 

Here we have defined an effective opto-mechanical driving amplitude VLm = go.-, and 



g = oj'zzp = uj' \J 2mlj optical cavity frequency shift per unit mechanical zero-point 

uncertainty. In the perturbative limit, this expression can be written in terms of the sus- 
ceptibility of a simple oscillator with an effective linewidth and frequency that is modified 
due to opto-mechanical interactions, x{^)~^ ~ ^mes — — i^jJ^eS- The effective linewidth 
and mechanical frequency shift can be interpreted as resulting from optically-induced cool- 
ing (or heating) and a dynamic optical spring constant, respectively. In the relevant regime 
of large detuning 6^u, n, and when the optical spring is dominant compared to the natural 
mechanical frequency, the effective mechanical frequency is given by 



UJm,eS^'2^m\l'^- (D3) 



The damping rate is given by 



tU-n 



1 



(«:/2)2 + {6 + u;^,eff)' {K/2y + {6- a;^,eff)^ 
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(D4) 



which is interpreted as the difference between anti-Stolces and Stokes scattering rates. Note 
that for positive detuning 6 > 0, the opto- mechanical interaction yields an increase in 
the mechanical frequency but an anti-damping force (r^s < 0). One can achieve simul- 
taneous stiffening and cooling by employing multiple beams with different amplitudes and 
detunings [H], but for our purposes it is sufficient to consider only the beam that leads to 
stiffening. 

We wish to consider the ratio of the effective mechanical frequency to the rate of decoher- 
ence induced by optical Raman scattering, which is given by the sum of the anti-Stokes 
and Stokes scattering rates. In the relevant limit of large detuning and dominant optical 
spring effect, one finds 

— — . p5 

This result states that the cavity must be driven very far off resonance in order to yield a 
frequency shift that is much larger than the decoherence rate. Operating at large detuning 
in turn requires extremely large cavity input powers to get an appreciable optical spring 
effect. As an example, we consider the dynamic spring constant for a realistic geometry, 
such as a Fabry-Perot cavity of length L = 1 cm and cavity finesse J-" = 10^ (with the cavity 
linewidth given by k = ttc/{J^L) = 27rxl50 kHz for our specific parameters). The optical 
driving amplitude is related to the input power Pi through fi^ = a/ kPi/ {2hwL) for perfect in- 
coupling efficiency, while the opto-mechanical interaction strength is of order u'^uq/L. The 
operating wavelength is taken to be A = 1 /im. We also assume that the SiN membrane has 
a radius of radius a = 10 fim and that it undergoes pure CM motion (such that its effective 
motional mass is the same as the physical mass). Then, an input power of Pj~2 kW is 
required if one wants to achieve a number of coherent oscillations N^"^^^ = ^^^J ^IO^ and an 
effective mechanical frequency of (X'm,cfT~27r x 1 MHz. This corresponds to an input intensity 
of ~10 W/yum^ for a beam focused to a size comparable to the membrane radius. 

In contrast, in our static trapping scheme, a comparable mechanical frequency and co- 
herence time can be achieved for an intra-cavity intensity of ~0.1 W/fim^, and the cavity 
can be driven resonantly to facilitate the intra-cavity field buildup. Because the static trap 
results in the membrane being trapped at an anti-node, there is no linear opto-mechanical 
coupling for the trapping field and the lowest-order opto-mechanical coupling is quadratic in 
nature. The decoherence rates caused by anti-Stokes and Stokes scattering in this case (at 
frequencies ujL±2um) have been calculated in Ref. [22] and are extremely rare for our real- 
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istic systems (occurring at a sub-Hz level). In fact, the dominant decoherence mechanism 
arising from the trapping beam is due to scattering out of the cavity and the photon recoil 
heating imparted on the membrane, as described later. 



Appendix E: Modification of opto-mechanical coupling strengths 

Our theory of optical trapping of membranes predicts dramatic corrections to the sim- 



ple model of opto-mechanical interactions given by Eq. (Dl), when the optical restoring 
forces (either static or dynamic) become large compared to the natural ridigity of the mem- 
brane. In this scenario of strong optical forces, the mechanical mode shape and thus the 
opto-mechanical coupling strength g become functions of intensity as well, with g generally 
decreasing with larger intensity. The origin of this effect is intuitively seen by considering 
a membrane that interacts with a Gaussian cavity mode whose beam waist w is smaller 
than the membrane radius a. Then, if the optical restoring forces are large compared to 
the membrane stiffness, the optical beam in fact resembles a new boundary condition that 
"pins" the region r < w of the membrane into place. This reduces the overlap between the 
mechanical displacement field and the optical beam, and thus g. 

This effect is illustrated in Fig. |4] for a free SiN membrane of thickness d = 30 nm and 
a = 25 /xm, interacting with a beam of waist w = 15 fim. In this calculation the membrane 
is statically trapped, although a similar effect would occur for sufficiently large dynamical 
backaction forces as well. In Fig. |4^, we calculate the opto-mechanical coupling strength 
g (to another cavity mode with the same beam waist but which exhibits an intensity gradient 
at the membrane position) as a function of the CM frequency, which is varied through the 
intensity of the trapping beam. The value of g is calculated using the expression [2j 

^oc z^pj dxdye-'^-"+y"^/'"\{x,y)/max\C\. (El) 

In Fig. |4^, we have normalized the obtained value of g with the value go if the motion were 
purely CM, where the displacement field ( is constant. At larger frequencies, g dramatically 
decreases, reflecting the "pinning" effect that the optical force has on the center of the 
membrane. This is also directly seen in Fig. lib, where we plot the displacement field C(x, y) 



for a CM frequency ujm/^ir = 300 kHz. Incidentally, Eq. (El) also makes it apparent that 



the torsional (1,0) membrane mode described in the main text has zero opto-mechanical 
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coupling, due to the odd versus even reflection symmetries of the torsional mode and cavity 
mode, respectively. 

Appendix F: Tethered membrane inside a Fabry-Perot cavity 

The membrane should be trapped within a Fabry-Perot cavity in order to enable optical 
cooling and reach the quantum regime [2D1 EI]- The membrane diffracts and scatters the 
cavity light, which introduces two important effects. First, the optical mode will no longer 
be Gaussian, and the new optical mode accommodated by the mirrors and the corresponding 
optical forces must be determined. Second, scattering out of the cavity reduces the cavity 
finesse, and the associated photon recoil acts as a stochastic force that heats the membrane 
motion. To calculate the cavity modes in the presence of the membrane, we use a modified 
Fox-Li propagation technique [23], as briefly described below. 

Specifically, the electric field is treated within the scalar paraxial approximation, and thus 
it is completely described by its transverse profile E{x,y). This approximation is justified 
by noting that the disk should primarily diffract light at small angles 0<{ka)~^ around the 
2;-axis, where k = 27t/X is the optical wavevector. Within this approximation, free propa- 
gation over a distance z is accounted for by a phase shift in the Fourier transform of the 
field profile, E{k^,ky)^e''''-'^^-+^i'^'/^'^^^E{k^,ky). In our case, we are interested in systems 
with rotational symmetry, and thus the transforms are implemented through the quasi- 
discrete Hankel transform described in Ref. |24J. Reflection off of a circular mirror with 
radius of curvature Rc and reflectance Rm is characterized by the real-space transforma- 
tion E{x^y)^\/R^E{x,y) ey.Y> (^ik{Rc— \JR1 — (x^ + Similarly, at the membrane 
location, the wave front experiences reflection and transmission amplitudes r{x,y), t{x,y), 
respectively, which multiply the incident real-space wave front there. For the case of a non- 
uniform disk, we approximate r{x,y),t{x,y) with the formulas for an infinite thin dielectric 
sheet of uniform thickness, replacing the uniform thickness d-^d{x, y) with the local thick- 
ness. Note that an initial wave front incident on the membrane thus splits into two wave 
fronts (a reflected and transmitted field), and we keep track of the multiple scattered fields 
to all orders in order to calculate the field buildup or cavity eigenmodes. In contrast, the 
original technique of Ref. [23j only accounts for transmission. Thus, our approach properly 
captures the effects of the reflected amplitude and back-scattered angle. Furthermore, our 
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modified technique reveals specifically at what frequencies resonances should occur. 

We now discuss the effect of the membrane on the cavity finesse. As realistic parameters, 
we consider a membrane placed symmetrically in the center of an optical cavity of length 
L = 1.99 cm with spherical mirrors having radii of curvature Rc = 1 cm and perfect 
reflectivity (such that the entire cavity linewidth k is attributable to scattering from the 
membrane). The transverse extent of the spherical mirror surfaces is = 0.95 mm, i.e., 
all portions of the beam front with + > are scattered out and set to zero upon 
reflection at the mirror. An empty cavity in this configuration yields a Gaussian mode of 
waist Wo~15 fim in the center. In Fig. [5^, we plot the membrane-limited cavity finesse 
Fmem = tic/kL for a membrane of uniform thickness d = 30 nm and varying radius a (black 
circles). Clearly, cavity losses are negligible when the nominal waist is small compared to 
the disk radius, wo/a<l. In the regime wo/a>l, however, the finesse rapidly drops, which 
is attributable to scattering by the hard edges of the disk. This effect is strongly reduced by 
"softening" or apodizing the disk edge [2S]. In Fig. [5^ (red circles), we also plot the finesse 
for a membrane whose thickness d{r) = do{l — {r/aYY tapers down to zero at the edge, 
where d^ = 30 nm is the maximum thickness. Remarkably, the apodization can improve the 
cavity finesse by several orders of magnitude. The modification of the cavity modes by the 
membrane is illustrated in Fig. ^p, where we plot the transverse profile at the membrane 
position for some representative apodized disk sizes. 



We find the CM eigenmodes of the apodized disk by using Eq. (A3), where the optical 
potential uJopt{r) is now evaluated using the modified cavity mode profiles. The thermoelastic 
limit is subsequently calculated using Eqs. (2) and (3) in the main text. In Fig. [sj:, the 
number of oscillations N^^^^^^ due to thermoelastic damping is plotted (in black). Here the 
circulating intra-cavity power is chosen such that the CM frequency is fixed at a;m/2vr = 
0.5 MHz. We next consider the effect of photon recoil heating. We assume that each 
scattered photon contributes the maximum possible momentum kick of hk along the z-axis, 
giving rise to a momentum diffusion process d{p1)/dt = {hk)'^Rsc [lOj, where Rsc is the 
photon scattering rate. Converting this expression into a jump rate, it can be shown that 
the number of coherent oscillations before a jump in the phonon number can be written as 

Here /max is the maximum cavity intensity (which in general does not need to be at the 
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center of the membrane due to mode distortion), coq = ck, V is the volume of the disk, 
and Vc is the cavity mode volume. Assuming that the cavity mode is not significantly 
distorted and that the beam waist WQ>a such that the entire membrane experiences the 

2 

optical force, one can approximate ^2^'^'"/^^ ^-*' Vc^ttWqL/A. In this case the number 
of coherent oscillations scales roughly as Note that this result is purely 

geometric in nature and also scales directly with the cavity finesse (which itself depends 
on V). In Fig. |5|3, we plot Nsc^^'' for the apodized disk (in red). Combining the effects 
of thermoelastic damping and recoil heating, the total number of coherent oscillations is 
given by N^Z^ = (iVir^"' + iV(°^'^-^)-i (blue curve). It can be seen that an apodized 
disk of radius r~9 /xm can support a coherence time of Nj:"^'^^ ^^2000 in a room temperature 
environment. 
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FIG. 1: a) Side view of a membrane supported by a single tether inside a Fabry-Perot cavity. The 
membrane has radius a and thickness d, while the tether has length L and a square cross-section of 
width b. It is trapped in the anti-node of a standing optical field with transverse intensity profile 
I{r) at the membrane location, b) Displacement fields of a few selected membrane modes (in 
arbitrary units), for zero trapping intensity. The black outline indicates the equilibrium position, 
(m, n) denote the number of nodal diameters and circles, respectively. The system dimensions are 
given by a = 10 ^m, b = d = 50 nm, and L = 50 fini. 
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FIG. 2: a) Normal mode frequencies of a free circular disk trapped in an optical standing plane 
wave, as a function of beam intensity. The disk has a thickness and radius of d = 50 nm and 
a = 10 /im, respectively, and material properties corresponding to stoichiometric silicon nitride, b) 
Normal mode frequencies of a rigid membrane suspended by a single tether, as a function of the 
optical restoring frequency Wopt acting on the membrane. The tether has length L = 50 /im and 
a square cross-section of & = 50 nm on each side, while the ratio of membrane to tether mass is 
given by M/rat = 125. Away from degeneracies, the mode spectrum consists of a CM mode with 
frequency ~a;opt and discrete tether modes with frequencies a;„ {n = 1, 2, 3, ...). Avoided crossings 
occur near degeneracies a;opt~'^n (see inset), c) Normal mode frequencies for a realistic tethered 
system, as a function of peak trapping beam intensity. The disk and tether have dimensions 
identical to those in a) and b), while the beam waist is it; = 35 /xm. The gray points indicate tether 
modes. The red (y)) green (' blue (□), and black (A) points denote the CM, (1,0), (2,0), and 
(0,1) membrane modes, respectively. 
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FIG. 3: Ratio f/opt/f^mech of optical to strain energy for the CM mode of different systems: a) 
an optically trapped free disk as a function of its frequency a;m/2'/r. Here the beam waist is fixed 
at t(; = 35 /im and the trap intensity is varied to yield the corresponding lo^- Inset: the trap 
intensity is fixed such that ct;m./2vr = 1 MHz in the plane wave limit while w/a is varied, b) a 
rigid membrane suspended by a single tether, as a function of the CM frequency, c) a realistic 
tethered structure, obtained by finite-element simulations. The system dimensions for these plots 
are identical to those in Fig. [2j namely, a = 10 /xm, L = 50 ^m, d = 6 = 50 nm. 




FIG. 4: a) Opto-mechanical coupling strength g oia, trapped free disk as a function of CM frequency. 
The coupling strength is normalized by the value corresponding to rigid (pure CM) motion qq. A 
decrease in g for increasing frequency is caused by the non-uniform optical force pinning the center 
of the disk in place. The dimensions for this simulation are d = 30 nm, a = 25 /im, and w = lb fim. 
b) Displacement field C{x, y) (in arbitrary units) for a free disk of the same dimensions, for a trap 
frequency of Um/'^'^ = 300 kHz. The displacement field clearly illustrates the pinning effect created 
by the optical forces. 
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FIG. 5: a) Membrane-limited finesse -Fmem of a Fabry-Perot cavity with a membrane in tfie middle. 
The black circles correspond to a flat membrane of uniform thickness d = 30 nm, while the red 
circles correspond to an apodized membrane with maximum thickness (Iq = 30 nm. The finesse is 
plotted as a function of the ratio of the empty-cavity beam waist wq to the membrane radius a. 
The cavity parameters are chosen such that wo = 15 ^m. b) Intensity profiles (in arbitrary units) 
of cavity mode in the presence of an apodized membrane. The intensity profile is evaluated halfway 
between the two cavity mirrors. The cavity and membrane parameters are provided in the main 
text. The blue and red curves correspond to disk radii a = wq and a = 2.5wo, respectively, while 
the dashed black curve is the Gaussian intensity profile for an empty cavity, c) The number of 
coherent oscillations of the CM motion of an apodized disk due to thermoelastic damping {N^^^^^ , 
black cirlces) and recoil heating {Nsc^'^\ red circles), as a function of disk radius and for fixed 
Wo = 15 fim. Also plotted is the total number of coherent oscillations {N^°^'^\ blue dashed curve), 
which is given by the sum in parallel of the individual contributions. 
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